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ABSTRACT 

The correlation between the soft X-ray and near infrared emission from AGN is 
analysed using composite models. We find new evidences for differences in ranges of 
parameters which characterize the NLR of Seyfert galaxies and LINERs. Soft X-rays 
show less variability, so they are better fitted for this kind of analysis. In our models 
soft X-ray are emitted in the post shock region of clouds with relatively high shock 
velocities Vs> 250 kms~^. Consequently dust emission peaks in the mid-IR. On the 
other hand, in the photoionized zone, dust is at lower temperature and usually does 
not contribute to the mid-IR emission. The results are sensible enough to allow the 
same modelling method to be applied to different types of AGN. We found that shock 
velocities are between 300 and 1000 kms^^, with the NLR of low-luminosity AGN 
and Seyfert 2 showing lower velocities than Seyfert 1 galaxies. The intensity of the 
ionizing radiation flux at the Lyman limit from the central source is low for LINERS 
and low-luminosity AGN (log Fh= 9 to 10), increasing towards Sy2 (log Fh'^ 11) 
and Syl (11 ^ Fh^l2). Results obtained by modelhng the Einstein and the ROSAT 
samples of galaxies are in full agreement. Dust-to-gas ratios by number are ^ 10^^'* 
in LINERs and LLAGN, between 10"^^ and 3 10"" in Syl and up to 5. lO^^^ in Sy2. 
In order to fit the infrared and X-ray continua, an rj factor is defined, which accounts 
for the emitting area of the cloud. If the infrared emission is due to bremsstrahlung 
and comes from the same cloud producing the soft X-rays, the 77 values obtained from 
both emissions must be the same. Therefore, if {vi)ir < (rf)aoft-x there must be a 
strong contribution of soft X-rays from the active centre. From the 77 values, we expect 
to identify the objects that could present strong variability. 

Key words: galaxies:Seyfert - radiation mcchanism:nonthermal - shock wave soft 
X-ray:galaxies - infrared: galaxies 



1 INTRODUCTION 

The main source of ionization and heating of the emission- 
line region in active galactic nuclei (AGN) is photoionization 
by radiation from the active centre (AC). Nevertheless, the 
contribution of shock waves, generated by the radial motion 
of the clouds or associated to jets, to the emission originated 
in the narrow line region (NLR) of Seyfert 1 (Syl), Seyfert 

2 (Sy2) galaxies, and LINERs was pointed out, for example, 
by Contini & Viegas (2001), Contini (1997), and Viegas & 
dal Pino (1992), and more recently by authors analysing the 
soft X-ray spectrum of AGN (e.g. Iwasawa et al. 2003). The 
relative importance of shocks and radiation can be roughly 
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evaluated in the different objects, because these processes 
act in different domains. In fact, shocks heat the gas within 
the clouds to high temperatures in the downstream region, 
so they can easily explain the observed high ionization lines 
(Rodriguez- Ardila et al. 2002) and soft X-ray emission (Con- 
tini, Prieto & Viegas 2002, and references therein). 

Indeed, if soft X-rays are emitted by bremsstrahlung 
from the NLR clouds, the gas temperature must be at least 
of the order lO" K. Such temperatures are obtained in the 
post-shock region if the shock velocity is higher than 250 
kms"\ The FWHM line profiles in the NLR of Seyfert 
2 and Seyfert 1 galaxies indicate that the cloud velocities 
can well be beyond this range. Modelling the emission-line 
and continuum spectra of some Seyfert 2 galaxies (e.g. NGC 
5252, Circinus, NGC 4151) indicate the presence of very high 
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shock velocities which can provide a thermal origin for the 
observed soft X-rays. 

Based on early observations of different galaxies by EX- 
OSAT, the nature and variability of soft X-rays was investi- 
gated, revealing a low and long term variability compatible 
with bremsstrahlung emission from a hot gas. For instance, 
Pounds et al. (1986) showed the existence of a soft X-ray 
spectral component, distinct from the main (nuclear) X-ray 
emission of NGC 4151 by a lack of variability. This compo- 
nent is shown to be consistent with thermal emission from a 
hot, confining medium in the region of the narrow emission- 
line clouds. Perhaps NGC 4151 is an exception to the general 
rule that soft X-ray emission from AGNs is variable, because 
most of its soft X-ray flux is absorbed by a high intrinsic 
column density, so that its faint, diffuse X-ray flux can be 
detected. This effect should not appear in most objects. 

The correlation between the soft X-ray and infrared lu- 
minosities in AGN has boon studied by several authors and is 
controversial. One important issue being the source of vari- 
ability in the various frequency ranges. In particular, the 
correlations of Lx (0.5-4.5 keV) versus L/h for the IRAS 
12/^m, 25/im, 60/im, and 100 /xm bands have been exam- 
ined by Green ct al. (1992). They claim that Lx versus 
L25^m and Lx versus Lioo/^m resemble very closely the 
Lx versus Leo^jm, correlation and that the plot of Lx versus 
L12 is considerably better defined, yet, just as in the Lso^im 
plot, "Seyfert galaxies still show no significant correlation" 
(Green et al., 1992, Fig. 3) 

On the other hand, the correlation between the 
10/xm and the hard X-ray emission showed opposite results 
(Carleton et al. 1987 and Giuiricin ct al. 1994) and was 
reanalysed by Giuricin, Mardirossian, & Mezzetti (1995). 
They obtained a wealc correlation with a harder (2 to 10 
keV) X-ray emission, mainly induced by Seyfert 2 data. 
They also found that the lO/xm emission correlates well with 
the 12/im and 25/im, and loss well with the far-infrared omis- 
sion at 60/im and lOO/^m. These authors favor thermal emis- 
sion as the origin of the mid-infrared emission of Seyfert 
galaxies. 

Dust is present in the NLR clouds and re-emission by 
grains may dominate the IR domain. If the intense UV- 
optical central radiation heating the nuclear dust has an 
X-ray counterpart, a correlation between the infrared and 
X-ray emission should be expected (Barvainis 1990). Be- 
cause in shocked zones the dust grains are usually heated to 
higher temperatures than in pliotoionizcd regions, the con- 
tribution from the shocked NLR clouds probably dominate 
the near infrared emission. A contribution to the IR from 
bremsstrahlung emission from a cooler gas is not excluded 
and would result in a correlation between the the soft X-ray 
and the near IR continua. Nevertheless, in the former case 
the mutual heating of dust and gas (Viegas & Contini 1994) 
can also lead to a correlation between those two continua. 
Therefore, in both cases a correlation is expected. Notice 
that the contribution of a dusty torus to the IR emission 
suggested by the AGN unified model (Antonucci & Miller 
1985) will give poor or no correlation with the soft X-ray 
emission as pointed out by Green et al.(1992) (see also Con- 
tini, Viegas, & Pricto 2003 and references therein). 

A major problem concerning correlations between X-ray 
and IR emission is X-ray variability. It is well established 
that soft X-rays from low luminosity AGN (LLAGN) and 



LINERs show less and/or longer variability than the hard 
X-rays (Halderson, Moran, & Filippenko 2001, Terashima, 
Ho & Ptak 2000). Following Green et al. (1992), we wiU, 
therefore, analyse the possible correlation between mid-IR 
(L12) and X-rays for these galaxies, which present a bolo- 
metric luminosity that can be several orders of magnitude 
lower than the brightest AGN. 

Notice that 10 /Ltm falls in the absorption trough of the 
spectral energy distribution due to silicate grains. On the 
other hand, PAH omission bands dominate the infrared SED 
of galaxies at 6.2, 7.7, 8.6, 11.3, and 12.7/xm (Spoon et al 
2002) . However, PAH grains are small compared with silicate 
grains and easily sputtered in AGN, where shock velocities 
arc high enough (> 300 kms"^). Therefore, we refer in this 
paper to the 12/itm emission in the mid-IR. 

It was found in previous investigations on the Seyfert 
galaxies, NGC 5252 (Contini, Prieto, & Viegas 1998a), Circi- 
nus (Contini, Prieto, & Viegas 1998b) and NGC 4151 (Con- 
tini, Viegas, & Prieto 2002), that the critical shock velocity 
(Vs) to produce dust emission in the mid-IR is of the order 
of 500 kms~^ (Viegas & Contini 1994), and that the crit- 
ical Vs for soft X-ray emission, fitting most of our sample 
for Seyfert 2 galaxies, is about 1000 km s~^ (Contini & Vie- 
gas 2000). However, the dust-to-gas ratio, d/g, may change 
from galaxy to galaxy, as well as in the different regions of 
a single galaxy. Therefore a scattering of the infrared data 
for a sample of galaxies is expected. 

The FWHM of the line profiles is lower in LINERs (< 
300 kms~^) than in Seyfert galaxies (300 to 500 kms~^). 
Notice, however, that broad Ha lines have been observed 
from LINERs (Ho, Fihppenko, & Sargent 1997), indicating 
that clouds with higher velocities may also contribute to the 
spectra. Therefore, in these objects, characterized by a low 
ionizing flux from the AC, shocks are relatively important 
and may contribute to both the near-IR and soft X-rays 
emission. The observed data will be confronted to model 
results obtained by rmmcrical simulations with the SUMA 
code (Viegas & Contini 1994, Contini & Viegas 2001), which 
accounts for both the shock and photoionization effects in 
the NLR clouds. 

Once tested for the low- luminosity objects, in the sec- 
ond part of this paper we will apply the same kind of anal- 
ysis to Seyfert galaxies. We will restrict the data to the soft 
X-ray part of the spectrum, where the variability is not as 
strong as in the hard X-rays. From the comparison with 
model results we expect to select those galaxies where the 
contribution to soft-X rays from the AC prevails, indicating 
in this way the objects of the sample which should show a 
higher degree of variability. 

The spectral energy distribution (SED) of dust and gas 
omissions obtained from the models, as well as the corre- 
sponding temperatures of the dust grains are discussed in 
Sect. 2. The Lx - L12 correlation for LLAGN and LINERs 
is investigated in Sect. 3, while that for Seyfert galaxies is 
discussed in Sect. 4. Concluding remarks appear in Sect. 5. 

2 GAS AND DUST EMISSIONS 
2.1 The model 

In our model the nuclear clouds are moving outward from 
the AC. A shock front appears at the cloud outer edge, while 



The Infrared - X-ray continuum correlation in AGN 3 




Table 1. Model results : the temperature of grains 



Figure 1. The SED of the continua emitted by bremsstrahlung 
and by dust from the clouds. Long-dashed lines refer to 
Vs=500 kms-i, logFh=12; short dashed lines to Vs=300 kms"!, 
logFh=ll. Solid lines refer to Vs=1000 kms~^. Thick lines rep- 
resent RD models, thin lines the SD models. The vertical lines 
correspond to the main IR observed wavelengths. 



radiation from the AC reaches the inner edge. The physical 
conditions of a emitting cloud due to the coupled effect of 
shock and photoionization by a radiation source arc obtained 
with the SUMA code assuming a plan-parallel symmetry. 
The calculations start at the shocked edge. In the first run, 
only the shock effect is accounted for, and the physical con- 
ditions are obtained for several 300) slabs before reaching 
the opposite edge. In this first run the results correspond to 
a shock-dominated model. Iterative runs alternatively start 
at the photoionized and shocked edges, accounting for both 
effects, until convergence of the physical conditions in the 
cloud is achieved. 

The input parameters for a single-cloud model are the 
shock velocity, Vg, the preshock density, no, the preshock 
magnetic field. Bo, the ionizing radiation spectrum, the 
chemical abundances, the dust-to-gas ratio by number, d/g, 
and the geometric thickness of the clouds, D. A power-law 
is adopted for the central ionizing radiation spectrum reach- 
ing the cloud, characterized by the flux Fh, at the Lyman 
limit (in units of cm~^s~^ eV~^) and the power law index 
a (FhOC For all the models, Bo = 10"'* gauss, auv 

= 1.5 and ctx ~ 0.4, and cosmic abundances (Allen 1973) 
are adopted. The model is either radiation dominated (RD) 
or shock dominated (SD) depending on the intensity of the 
ionizing flux and on the shock velocity (Viegas-Aldrovandi 
& Contini 1989). For SD models Fh= 0. 

The gas entering the shock front is thermalized and 
heated. Depending on the shock velocity, a high tempera^ 
ture zone can be produced where soft X-rays are generated 
by bremsstrahlung. A lower temperature zone may produce 
near infrared emission by the same mechanism. The tem- 
perature of the photoionized zone is never too high (< 4. x 
10* K), thus its bremsstrahlung emission is mainly in the 
IR-optical range, depending on the intensity of the ioniz- 
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ing radiation. Thermal emission from dust grains could also 
produce near IR emission, as long as the dust temperature 
is high enough as discussed below. 

Before discussing the temperature of the grains, let us 
illustrate the SED of the continuum accounting for both gas 
bremsstrahlung and dust re-emission calculated by a com- 
posite model (Contini & Viegas 2000). The SED shown in 
Fig. 1 correspond to single-cloud models with various cloud 
velocities and ionizing radiation intensities, assuming a ge- 
ometric thickness of the cloud D = 10^^ cm and dust-to-gas 
ratio d/g =10"^'' by number (or 4. lO""* by mass, adopting 
silicate grains). Following the results obtained in our previ- 
ous papers, models with high velocities (1000 kms~^) are 
SD, because the effect of the shock prevails on photoioniza- 
tion from the AC. The results shown in Fig. 1 will help to 
understand the soft X-ray emission radiation from clouds in 
the physical conditions suitable to those in the NLR of AGN 
and in LINERs. 



2.2 The temperature of the dust grains 

The dust temperatures calculated by models with Vs= 300 
and 1000 kms"'^ are presented in Table 1, for both SD side 
and the RD side of the the single-clouds. We assume that 
grains are present inside the whole cloud before entering the 
shock front and all the dust grains have initially the same 
The grain radius may change due to sputtering. A 
dust-to-gas ratio d/g = 10~^* by number is adopted. 

In the photoionized zone, sputtering is not important 
and the grain temperature increases with decreasing grain 
radius and increasing Fh. The behavior is rather different 
in the shocked zone because of sputtering throughout the 
shock front. The higher Vs the stronger the sputtering and 
the higher the grain temperature, as long as the grains sur- 
vive sputtering. So, larger grains, which survive sputtering 
through high velocity shock fronts, show a distribution of 
different radius throughout the cloud. They reach higher 
temperatures and the peak of dust thermal emission shifts 
toward higher frequencies. This peak can be high depending 
on the volume of the dusty region. When the grain initial 
radius is small and the shock velocity high the dust grains 
are destroyed by sputtering, so dust cannot survive in the 
SD side of the cloud (see Contini, Viegas, & Pricto 2003). 
Moreover there is a mutual heating between gas and dust 
(see Viegas & Contini 1994). 

It can be seen that Tg,. in the RD zone is always lower 
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Figure 2. The ratio of tlie 12 fj, luminosity to the soft X-ray 
luminosity versus the soft X-ray luminosity for the low-luminosity 
active galaxies of the Halderson et al (2001) sample. 

than in the SD zone, even for small grains (agr=0.01 ^m) 
and high Fh (5 10^^ cm"^s~^eV~^). On the other hand, 
dust is destroyed by sputtering if the initial grain radius is 
low and Vsis high, as exemplified by the model with a,gr= 
0.01 /xm Vs= 1000 kms"^ 

In brief, if the mid-infrared emission of AGN, particu- 
larly the 12jj, emission, is thermal emission due to grains, 
the grains must be heated in a shocked zone. 



3 LOW LUMINOSITY AGN 

The correlation between the near IR and the soft X-rays 
for LLAGN is analysed using 36 nearby objects from the 
Halderson et al. (2001) sample. Most of the X-rays data 
come from observations with ROSAT while the infrared data 
from IRAS. Some of the objects are also included in the 
Terashima ct al. (2002) sample of LLAGN observed with 
ASCA. Short term variability ( ^ 1 day) was found in NGC 
3031, NGC 4258 and NGC 5033. Long-term variability is 
observed in NGC 3031, NGC 4258 and NGC 4579, while no 
significant variability was found in NGC 404, NGC 4111, 
NGC 4438, NGC 4579 and NGC 4639. 

The correlation between the 12 /xm and soft X-rays for 
the LLAGN is shown in Fig. 2, where we plot the ratio of the 
IR and X-rays observed fluxes versus the X-ray luminosity. 
The choice of the vertical axis is explained below, where we 
compare the observational data with the results of numerical 
simulations. 

3.1 Modelling the sample 

The single-cloud model provides the continuum and 
emission-line fluxes calculated at the cloud. Thus, in order 
to obtain the luminosities, the model results must be mul- 
tiplied by a factor 77 accounting for the emitting area. As- 
suming that the single cloud conditions correspond to the 



Figure 3. The comparison of model calculations with the obser- 
vational data. The results of the models with different values of 
Vsand the same log Fjjajre joined by dotted (log Fij= 9), short- 
dashed log Fh= 10), long-dashed ( logFh=ll) and dash-dotted 
(logFij= 12) lines. Solid lines correspond to shock dominated (SD) 
models calculated with Fij= Thin lines refer to bremsstralilung 
emission at 12 /im, while thick lines refer to dust emission at 12 
//m. The values of Vsincreases from left to right. 

gas around the central source, the emitting area can be re- 
lated to the distance of the cloud to the central source and 
to the covering factor (see, for instance, Contini & Viegas 
2000). Plotting a flux ratio versus a luminosity, as in Fig. 2, 
this factor will affect only the horizontal axis, and the match 
of the models to the data would be easier, if IR and X-ray 
emission come from the same cloud and there is a negligible 
contribution of X-rays from the central source. 

Since we assume that the soft X-rays come from the 
shocked edge of the NLR clouds, the post-shocked zone 
must reach temperatures higher than 10^ K and a mini- 
mum shock velocity of 300 kms"'^ is required. Based on 
previous works, a grid of models with Vs=300 kms~^and 
no=300 cm-3, Vs=500 kms'^and no=300 cm'^, Vs=700 
kms~^and no=500 cm~^, and log Fh= 9,10,11,12 was ob- 
tained for clouds with a dust-to-gas ratio of 10^^^, a geo- 
metrical width of D=10^^ cm, a pre-shock magnetic field of 
Bo=10~^ gauss and initial grain radius of 0.2 /xm. Notice 
that log Fh= 12 is an upper limit since a low radiation flux 
from the AC is characteristic of LLAGN. 

Model results and observational data are compared in 
Figure 3. The lines join results of models with different val- 
ues of Vsbut the same Fh. 

The observations of the galaxies NGC 4486, NGC 4579, 
and NGC 4051 are better explained by SD models (thin 
solid line joining the results for Vs= 900, 1000, and 1100 
kms~^ and d/g = lO"'^^). Only the bremsstrahlung is shown 
in this case, because at these high velocities the grains 
are mostly destroyed by sputtering. The Fi2;jm /Fx ra- 
tios for SD models calculated with Vs= 300, 500, and 700 
kms~^ and d/g=10~^* are too low. Higher values (thick 
solid line) are obtained assuming d/g = 2xl0~^^. 
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Table 2. The characteristics of LLAGN 
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The observed correlation is reproduced by models cal- 
culated with log Fh=9 and 10. However, it is difficult to 
distinguish between bremsstrahlung and thermal emission 
by dust. 

Assuming all clouds at the same distance (r) from the 
center with the covering factor equal to unity, the best gen- 
eral fit to the observations, which reproduces the slope of 
the correlation, is obtained with clouds at r = 15.45 pc. 



3.2 Results for individual galaxies 

Single-cloud models for the galaxies in the Halderson et al. 
sample are chosen based on the two data: the infrared lumi- 
nosity at 12 /xm and that in the soft X-rays (0.5-2.4 keV). In 
the latest years point-to-point observations of a single ob- 
jects permit to obtain information about the distribution of 
velocities, densities, intensity of the radiation from the AC, 
and dust-to-gas ratios throughout the galaxy. By modelling 
the soft X-ray continuum, however, we select the clouds with 
shock velocities high enough to heat the gas in the post- 
shock region to high temperatures. The modelling of these 
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Figure 4. The best fit of the data by model results (see text). 



clouds permits to evaluate roughly their distances from the 
AC. In Table 2, for each galaxy of the sample the average 
characteristics of the emitting clouds reproducing the ob- 
served correlation are listed. 

The galaxies of the Halderson et al. sample are sepa- 
rated in LINERs (empty squares) and LLAGN (full squares) 
in Fig. 4. Some of the curves representing the results from 
models shown in Fig. 3 appear horizontally shifted in Fig. 
4 in order to reach most of the objects. This shift corre- 
sponds to an T) factor different from that used in Fig. 3. 
As seen in Fig. 4 most of the LINERs are fitted by low 
Fh (log Fh=9-10), confirming that low Fh are a characteris- 
tic of these objects (Contini 1997). The clouds in NGC 404 
and NGC 4395 that arc closer to the AC, arc reached by a 
jather high flux (log Fh=ll). Our results suggest that NGC 
4168 should be classified as a LINER. NGC 3031 has shock 
dominated clouds close to the nucleus (r=2.4 pc) heated by 
high velocity shock (Vs=700 kms-^). NGC 1068, on the 
other hand, shows IR emission by dust from clouds with 
d/g=10-i* and Vs=700 kms-\ located at 23.6 pc from the 
center. Clouds in NGC 4486 and NGC 4579 have the highest 
shock velocities (Vs~ 1000 kms-^) and are located in the 
nuclear region (r=8.9 pc). Notice that NGC 1068, which has 
an Ha luminosity greater than the limit for LLAGN (10*" 
ergs s-i), does not follow the general trend (see Fig.2), since 
high velocity clouds at large distance from the center are re- 
quired to account for the observed Ljc (Table 2). NGC 4051 
is a warm absorber and soft X-ray emission has been ex- 
plained by a more complete model (Contini & Viegas 1999). 



3.3 Composite radio emission 

Before concluding the analysis of the LLAGN, we would like 
to strengthen the fact that shocks and the ionizing radiation 
from the AC act together heating and ionizing the gas also 
in low luminosity objects as illustrated by the radio con- 
tinuum. In fact, the trend of the SED in the radio domain 
is obtained by summing up the contribution from the syn- 
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Figure 5. The calculated radio continuum for NGC 4374, 
NGC 4278, and NGC 4258. The two components are shown: 
the power-law due to Fermi mechanism (thin solid line) and the 
bremsstrahlung from a cool gas (dotted line) The thick solid line 
corresponds to the sum of the two components. 



chrotron radiation created by the Fermi mechanism at the 
shock front and the bremsstrahlung from the cooler cloud 
zones. 

A few examples appear in Fig. 5, showing model re- 
sults and observational data for the galaxies NGC 4374 
(top panel, NGC 4278 (middle panel), and NGC 4258 (bot- 
tom panel) . The model input parameters are listed in Table 
2. The observational data were found in the NASA Extra- 
galactic Data (NED)) and come from Laing & Peakock 1980, 
Kuhr ct al. 1981, Kollcrmann ct al 1969, Gowcr ct al. 1967, 
Wright & Otrupcek 1990), Douglas et al. 1996, Ekers 1969, 
Large et al. 1981, Heeschen & Wade 1964, White & Becker 
1992, Waldram et al. 1996, Becker ct al. 1995, Condon ct al. 
1983, Dressel & Condon, 1978, Becker et al. 1991, Gregory 
& Condon 1991, Israel & Mahoney 1990, Ficarra et al. 1985. 



4 THE Lx - L12 CORRELATION FOR 
SEYFERT GALAXIES 

In this section the correlation between the soft X-ray lumi- 
nosity {Lx ) and the 12 /um luminosity (L12) is investigated 
for a relatively large sample of Scyfcrt galaxies which in- 
cludes objects with higher luminosities relative to those of 
LINERS and LLAGNs studied in the previous sections. 

We have collected a sample of 31 Syl and 19 Sy2 ob- 
served with the Einstein (0.5 - 4 kcV band) and IRAS satel- 
lites. Another sample includes 26 Syl, 8 Syl. 5 and 10 Sy2 
observed with ROSAT (0.5 - 2.4 keV). The Einstein data 
come from Green, Anderson & Ward (1992), the ROSAT 
data were taken from Moran, Halpern & Helfand (1996) . The 
IRAS data come from the NED. The luminosities arc calcu- 
lated assuming a Hubble constant equal to 50 kms^^Mpc"^. 

The correlation between L12/Lx and Lx for Syl and 
Sy2 galaxies are shown in Fig. 6 using Einstein (top panel) 
and ROSAT data (bottom panel), respectively. For com- 
parison, the LLAGN data (Helderson et al. 2001) are also 
plotted. As expected the LLAGn are fainter in X-rays than 
the Seyfert galaxies. As expected, NGC 1068, included in 
the Halderson et al. sample, is located in the Seyfert galaxy 
zone. 

The separation between the two classes of objects is 
more clearly marked for ROSAT data. The L12/Lx ratios 
are roughly of the same order. These ratios depend mainly 
on the dust-to-gas ratio and on the physical conditions in 
the NLR, indicating that the conditions are similar in the 
two classes of galaxies. 

The same series of models (same notations) applied to 
the LLAGN in Fig. 3 are used to roughly model the Seyfert 
galaxies. In order to fit the data for the Einstein sample, we 
adopted an r; corresponding to an average distance of the 
clouds from the AC of ~ 460 pc for Syl and ~ 90 pc for 
Sy2. This corresponds to shifting the model ensemble to the 
right hand side in the diagram. However, to find the best 
fit for Sy2 wc had to shift the models upwards by a factor 
of 30. The vertical-axis corresponds to the L12 / Lx ratio, 
where L12 represents dust emission. Therefore, the best fit 
corresponds to dust-to-gas ratios higher by a factor of 30 
(d/g = 3 X 10"^'^) and already suggests that the NLR of Sy2 
is dustier than that of Syl. The same parameters used for 
the Einstein sample are adopted to fit the ROSAT sample 
(Fig. 6 bottom). Notice that Mrk 335, NGC 1365, and NGC 
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Table 3. The logarithm of the rj factors using the Einstein data 



type model X-ray^ IR(br.) IR(dust) r (pc) d/g (IQ-i*) line type 



Syl 


RD 


43.8 


44.2 




747. 




thick long-dashcd 


Syl 


SD 


42.6 




46.0 


187.6 


251. 


thin short-dashed 


Syl 


RD 


44.1 




45.2 


1055. 


1.26 


thick short-dashed 


Sy2 


SD 


40.1 




44.5 


10.5 


2512. 


thin short-dashed 


Sy2 


SD 


40.1 




45.3 


10.5 


1.58(4) 


thin dotted 


Sy2 


RD 


43.3 




44.8 


420. 


3.16 


thick short-dashed 


Sy2 


RD 


43.3 




45.3 


420. 


10.0 


thick dotted 



Einstein data 



Table 4. The logarithm of the rj factors using the ROSAT data 



type model X-ray^ IR(br.) IR(dust) r (pc) d/g (10"") line type 



Syl 


RD 


44.1 


44.1 




1055. 




thick long-dashed 


Syl 


SD 


42.5 




46.2 


167.2 


501. 


thin short-dashed 


Syl 


RD 


43.9 




44.9 


838. 


1. 


thick short-dashed 


Sy2 


SD 


40.1 




44.3 


10.5 


1585. 


thin short-dashed 


Sy2 


SD 


40.1 




45.1 


10.5 


1.(4) 


thin dotted 


Sy2 


RD 


43.2 




44.6 


374. 


2.5 


thick short-dashed 


Sy2 


RD 


43.3 




45.4 


420. 


12.6 


thick dotted 



1 ROSAT data 



1566 from the Einstein sample, showing low Lx , should be 
classified as Sy2 rather than Syl. 

The fit shown in Fig. 6 is reasonable, although more 
refined modelling is necessary. This will be done in the next 
section, where the NLR physical conditions, as well as the 
evaluation of the distances of the Syl and Sy2 emitting 
clouds to the center, are discussed. For this, the L12 ver- 
sus Lx correlation is used for both Syl and Sy2 galaxies 
(Fig. 7 and 8, respectively, corresponding to the Einstein 
and ROSAT data). 

4.1 Modelling the correlation 

Single-cloud models that better reproduce the L/ji - Lx cor- 
relation shown in Figs. 7 and 8 are chosen from a grid of 
models. All models are calculated with d/g — 10~^^ by num- 
ber, Bo= 10"* gauss, and no = 300, 500, and 1000 cm~^ for 
models with Vs= 300, 500, and 1000 kms~^, respectively 
(cf Fig. 1). Usually in the NLR the decreasing radiation flux 
with the distance to the AC is followed by a decreasing ve- 
locity field. Moreover, wc adopted a low dust-to-gas ratio as 
a lower limit, in order to find out the whole range of d/g 
values for both Syl and Sy2. 

The theoretical results are represented by lines, labelled 
by the input parameters Vs for SD models (thin lines), and 
the pair (Vs,log Fh) in the case of RD models (thick lines). 
Short dashed Unes refer to dust emission in the IR, while 
long dashed lines to IR bremsstrahlung. 

Fitting of the observed correlation with our models is 



not so simple. First, recall that models correspond to single 
clouds (sec Sect. 3.2). Thus, the comparison between the 
theoretical and observational data can only give a rough 
indication of the physical conditions of the gas and of the 
velocity field of the clouds emitting soft X-rays. Second, as 
already commented above (Sect. 3.1), in order to derive the 
luminosities, the fluxes calculated by the models must be 
multiplied by a factor (rj) which represents the fraction of 
the observed luminosity coming from a specific type of cloud, 
characterized by Vs, no and Fh. If the infrared emission is 
due to dust, this factor must also account for the differences 
of the dust-to-gas ratio assumed in the calculations and the 
real value for each galactic nucleus. A higher d/g value leads 
to an enhancement of the L/_r emission. Thus, we expect 
that the 77 factor is not constant but may have different 
values for different galaxies. It could be empirically obtained 
by fitting the observed emission-line and continuum spectra 
with a multi-cloud model, as already presented for Circinus 
and NGC 5252 (Contini et al. 1988 b,a). 

Since our goal is to analyse the L12 - Lx correlation we 
try to obtain a rough picture of the physical processes in the 
active nuclei with a minimum number of input parameters. 
First we constrain the models by considering the main trends 
shown by the data, because these trends define which type 
of models must be adopted (RD , SD, IR from dust or IR 
from bremsstrahlung). Then, we shift the model results in 
order to cover the largest number of data by changing the r; 
factor. 

Recall that in the SED presented in Fig. 1, the high 
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Figure 6. Top panel : L12/Lx versus Lx for Syl (filled cir- 
cles) and Sy2 (filled triangles) galaxies using Einstein data. Bot- 
tom panel : L12/Lx versus Lx for Syl (open circles) and Sy2 
(open triangles) galaxies using ROSAT data. In both panels, filled 
squares represent the LLAGN and the lines correspond to model 
results with the same notation used in Fig. 3. 



frequency domain corresponds to the hot gas emission, i.e., 
even in radiation dominated models it depends on the shock 
velocity. So, the soft X-ray emission is due to bremsstrahlung 
in the shocked zone of high velocity clouds while its infrared 
counterpart can contribute to the mid- or near-IR emission. 
In this case, since both continua come from the same cloud, 
the corresponding r; should be the same. 

The logarithm of r? for the bremsstrahlung and for dust 
emission, adopted for modelling the correlations presented 
in Figs. 7 and 8, arc given in Tables 3 and 4 for Einstein and 
ROSAT data, respectively. In the first column of Tables 3 
and 4 the type of the galaxy is given. In the ROSAT sample, 
Syl includes the Syl. 5 data. The model type (SD or RD) is 



indicated in the second column. The logarithm of the r) value 
which provides the best fit of the data by the corresponding 
model (see line type in column 8) are given in columns 3, 
4, and 5 for Lx , Ij12 generated by bremsstrahlung, and 
L12 by dust, respectively. In columns 6 and 7 the calculated 
distance of the emitting clouds from the AC (assuming a 
covering factor equal to unity) and the d/g ratio by number 
are given, respectively. The values of d/g by number which 
better fit the data are calculated from the difference between 
r/ corresponding to IR(dust) and r/ corresponding to X-ray 
for each type of galaxies in Tables 3 and 4. The line type 
used to draw the curves corresponding to the models (Figs. 
7 and 8) are listed in the last column. 



4.2 Seyfert 1 galaxies 

Two clear trends are present in Fig. 7 (top panel) corre- 
sponding to the Syl Einstein sample: (a) Syl with Lx ^ 3 
X 10*^ show a larger scatter of the L12 values; (b) the ob- 
jects with higher Lx values are mainly located between RD 
(thick short-dashed line) and SD ( thin short-dashed line) 
models both with IR due to dust emission . The results cor- 
responding to L12 due to bremsstrahlung (thick long-dashed 
line) could eventually explain some Syl and Sy2 data. 

Regarding the ROSAT data (Fig. 7, bottom panel), 
most of the data arc also inside a sector defined by RD mod- 
els (the low edge) and SD models (upper edge), as found 
for the Einstein sample. High L12 Syl galaxies are bet- 
ter fitted by dust emission from SD models with 300 < 
Va^ 1000 kms~^ (thin line), whereas those with low L12 
would be dominated by dusty RD clouds with Fh= 10^^^^^ 
( cm-^s"^eV"^) and Vs in the range 300 - 500 kms"^ 
Also for the ROSAT data, the line corresponding to the IR 
bremsstrahlung emission fits some Syl and Sy2 data, but 
not the trend of the Syl complex. 



4.3 Seyfert 2 galELxies 

The sample of Sy2 is shown in Fig. 8 both for Einstein (top 
panel) and ROSAT data (bottom panel). In both samples, 
the data are better explained by RD or SD models with 
L12 produced by dust thermal emission. The notation is the 
same as in Fig. 7. Notice that thin short-dcish lines refer 
to SD models and thick short-dash to RD models with IR 
from dust, while dotted lines correspond to the same mod- 
els shifted in the diagrams toward higher L/ij , in order to 
explain the high IR luminosity of objects which are rich in 
dust. 

Except for four objects in the top panel and three in the 

bottom panel, most of the data points are located between 
two SD curves, characterized by log r; = 40.1 and 43.3 On 
the other hand, the objects on right of the the SD curves are 
close to the RD results (thick short-dash) with log Fh=ll, 
Vs between 300 and 500, and a relatively low d/g (see Tables 
3 and 4). 



4.4 Discussion 

The first results of our investigation show that the X-ray 
and near-IR correlation for Syl can be explained by the 
composite models of clouds producing both bremsstrahlung 
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Figure 7. L12 versus Lx for Seyfert 1 galaxies (larger symbols). 
Top panel refers to the Einstein sample, while the bottom panel 
to ROSAT sample. The curves are labelled by the Vsand log 
Fhvalues for RD models and by Vsfor SD models. The line types 
are given in Tables 3 and 4 (last column). For the Einstein data 
the notation is: Sy 1 (filled circles) ; For the ROSAT data the open 
circles represent Syl, and small open circles the Sy 1.5. In both 
diagrams, tiny symbols correspond to Sy2 data. 



Figure 8. L12 versus Lx for Seyfert 2 galaxies (larger sym- 
bols). Top panel refer to the Einstein sample, while the bottom 
panel to ROSAT sample. The curves are labelled by the Vsand 
log Fh values for the RD models and by Vsfor SD models. The 
line types are given in Tables 3 and 4 (last column). Filled trian- 
gles represent the Einstein sample (top panel) and open triangles 
(bottom panel) the ROSAT sample. In both panels tiny symbols 
correspond to the Syl data. 



and dust emission. In the diagrams these models cover most 
of the area where the AGN data points are distributed. 

In a photoionized zone, the dust cannot reach a tem- 
perature high enough to contribute to the 12/im emission, 
whereas a shocked zone may heat the dust to higher temper- 
atures, as shown in Table 1. Therefore, models accounting 
for the shock are important. Considering the whole sample 
for a given class (Syl or Sy2), the contribution of RD and 
SD clouds to the IR emission coming from thermal emis- 
sion by dust provide the best fit to the correlation for both 
Einstein and ROSAT data. 

The values of the factor 77 adopted for the models and 



presented in Tables 3 and 4 give a rough estimation of the 
distance of the emitting nebula from the AC. Notice, how- 
ever, that the observed Lx may include a contribution from 
the central source, while L12 depends on d/g. So, the 1) 
value corresponding to L12 due to bremsstrahlung, which is 
related to the distance of the cloud to the AC, can be lower 
than that corresponding to Lx • 

The ratio of r\ for Lx and L12 will provide an estimate 
of the Lx emitted from the AC relative to Lx emitted by 
bremsstrahlung. For example, the observational data (Ein- 
stein sample) for the Syl galaxy Mrk 335, could be repro- 
duced by IR from dust with a high d/g. However, it could 
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be also reproduced by IR produced by bremsstrahlung, with 
\og{'ri)soft-x='i3.7 and \og{ri)iR—43.6. A difference in the r] 
values showing \og{ri)soft-x higher than log(?7)iH indicates 
a contribution from the AC to Lx • In this case, this galaxy 
could show soft X-ray short term variability. 

The theoretical results shown in Figs. 7 and 8 corre- 
spond to models with d/g =10~^^. With a lower d/g, L12 
dust will haxdly prevail over L12 bremsstrahlung, therefore, 
also in the case of emission by dust we assume as indicative 
of the nebula radius the lowest value of rj. The Sy2 galaxies 
included between the thin short-dash and dotted curves of 
Fig. 8 have a d/g between 2.5 10"^^ and 1.6 10"" (Einstein) 
and between 1.6 10"^^ and 10"" (ROSAT). The results for 
the Syl galaxies indicate d/g = 2.5 10"^^ - 1.3 10"" (Ein- 
stein) and 5. 10"" - 10"" (ROSAT). 

Regarding the objects with IR coming from dust re- 
radiation, results obtained for the galaxies of both Einstein 
and ROSAT samples, we obtain from the values of rj listed 
in Tables 3 and 4 a characteristic distance of a few hundred 
pc to a few Kpc for the Syl. The distances for Sy2 galaxies 
are smaller going from a few pc to about 500 pc. 

In brief, the results obtained by the modelling the Lx - 
L12 correlation suggest that both SD and RD clouds arc 
present in Syl and Sy2. Notice, however, that for Sy2 galax- 
ies, Va is about 300 km s"^ and log Fh is of the order of 11, 
while Syl have higher Vg (up to 600 kms~^) and higher 
Fh (log Fhup to 12). 

Clouds closer to the nuclear region are shock dominated 
because a large amount of dust can prevent the ionizing ra- 
diation from the AC. High velocity (Vb=1000 kms~^) and 
dust rich SD clouds are more likely to be present in Sy2 than 
in Syl, and are closer to the center. This indicates that the 
emitting region in Sy2 is more compact than in Syl. How- 
ever, the maximum distance of a few Kpc obtained for Syl 
from the fit of the ROSAT sample could represent an upper 
limit if the omitting gas has both a high d/g and a strong 
contribution to the X-ray coming from the active nucleus. 
Schmitt et al. (2003) found that both Seyfert types have 
similar distributions of the NLR sizes, although it is ex- 
pected that Syl have in average smaller NLR sizes than Sy2 
(Mulchaey ot al 1996). Notice that Schmitt ct al claim that 
there are "effects which could alter both the measurements 
and models" . 

Finally, the results obtained from the ensemble of the 
Seyfert galaxies reveal that the X-ray - near-IR emitting 

regions of Syl and Sy2 have indeed different characteristics, 
which are hardly noticeable when modelling single objects. 



5 CONCLUDING REMARKS 

In this paper we looked for new evidences for differences in 
parameter ranges characterizing the NLR of Syl, Sy2, and 
LINERs (and LLAGN) by modelling of the Lx - L12 corre- 
lation by composite models accounting for dust emission. 

Soft X-rays are considered because they show less dra- 
matic variability. In our models, soft X-ray are emitted in 
the post-shock region of clouds with relatively high shock ve- 
locities. Dust and gas are coupled crossing the shock front, 
therefore, dust emission from the cloud region that emits 
the soft X-rays peaks in the mid-IR. Shocks are relatively 
important in modelling the LLAGN spectra, so, we have 



started our investigation with the L/_r - Lx correlation for 
LLAGN. The results are sensitive enough to apply the same 
modelling method to Seyfert galaxies. 

We found that shock velocities are between 300 and 
lOOO kms~^for all the objects, lower in LLAGN and Seyfert 
2, higher in the NLR of Seyfert 1 galaxies. The intensity of 
the flux from the AC is low for LINERs and LLAGN (log 
Fh=9 - 10) and increases towards Sy2 (log Fh— 11), being 
higher for Syl (log Fh<12). Results obtained by modelling 
the Einstein sample and the ROSAT sample axe in full agree- 
ment. Dust-to-gas ratios by number axe laxger than 10"^'' 
(4. lO"-* by mass) for LINERs and LLAGN, between 10"" 
and 3 x 10"^'^ (4. 10""* and 0.012 by mass) for Syl, and 
reach 2. x 10"" (0.8 by mass) for Sy2. 

The soft X-ray variability of Seyfert galaxies can be 
inferred from the r] values deduced from the models. If {ri)iR 
< {'n)soft-x there is a strong contribution of soft X-rays 
from the AC to the observed value. Therefore, objects of 
the sample used in this paper which could present strong 
variability can be identified. 
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